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Abstract Although rare, collisions of two or more bodies in the A^-body problem are 
apparent obstacles at which Newton's Law of Gravity ceases to make sense. With- 
out understanding the nature of collisions, a complete understanding of the N-body 
problem can not be achieved. Historically, several methods have been developed 
to probe the nature of collisions in the A^-body problem. One of these methods re- 
moves, or regularizes, certain types of collisions entirely, thereby relating not only 
analytically, but also numerically, the dynamics of such collision motions with their 
near collision motions. By understanding the dynamics of collision motions in the 
regularized setting, a better understanding of the dynamics of near-collision motions 
is achieved. 



1 Introduction 

For ages, humankind has observed the regular and predicable motion of the planets 
and other bodies in the solar system and asked, will the motion of the bodies in the 
solar system continue for ever as they are currently observed? This philosophical 
question is the object of the mathematical notion of stability. A difficulty in applying 
the notion of stability to the motion of the solar system is that of collision and 
near-collision motions of bodies in the solar system. Collision and near-collision 
motions do occur in the solar system. Section [2] recounts a few of these that have 
been observed or predicted. 

The standard mathematical model for understanding the motion of planets and 
other bodies in the solar system is the Newtonian A^-Body problem, presented in 
Section [3] In included here are some of the basic features and mathematical theory 
of the Newtonian A^-Body Problem, its integrals or constants of motion, special 
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solutions such as periodic solutions, and the notions of stability and linear stability 
of periodic solutions and their relationship. 

The notions and basic theory of collisions and singularities in the Newtonian 
Af-Body Problem is presented in Section [4] This includes a discussion of the prob- 
abilities of collisions, and the regularization or the lack thereof for collisions. A 
collision motion is rare in that is has a probability of zero of occurring, whereas 
a near collision motion has a positive probability of occurring. Regularization is a 
mathematical technique that removes the collision singularities from the Newtonian 
Af-Body Problem and enables an analysis of near-collision motions in terms of col- 
lision motions through the continuous dependence of motions on initial conditions. 
This regularization is illustrated in the collinear 2-Body Problem, the simplest of all 
the A^-Body Problems. 

Recent results are presented in Section[5]on the analytic and numerical existence 
and numerical stability and linear stability of periodic motions with regularizable 
collisions in various A^-Body Problems with N = 3 and N = A. Although fictitious, 
these periodic motions with regularizable collisions provide a view of their near- 
collision motions which could be motions of the bodies in the A^-Body Problem that 
are collision-free and bounded for all time. 



2 Phenomenon 

Collisions and near-collisions of two or more solar system bodies are apparent ob- 
stacles at which Newton's Law of Gravity becomes problematic. Velocities of col- 
liding bodies become infinite at the moment of collision, while velocities of near- 
colliding bodies become very large as they pass by each other. Both of these situa- 
tions present problems for numerical estimates of the motion of such bodies. 

Although collisions are rare, historical evidence of collisions of solar system 
bodies is viewable on the surface of the Earth and the Moon [8 |. Only recently have 
collisions of solar systems bodies actually been observed. As the comet Shoemaker- 
Levy 9 approached Jupiter it was torn apart into fragments by tidal forces. In July of 
1994, at least 21 discernible fragments of Shoemaker- Levy 9 collided with Jupiter. 
These were the first ever observed collisions of solar system bodies. An animation 
of some of the fragments of Shoemaker-Levy 9 colliding with Jupiter can be found 
at w w w2 jpl .nasa. gov/sl9/anim. html. 

Near-collision motion are less rare than collisions. As of March 2012, there are 
nearly 9000 known near-Earth asteroid^ of which 1306 are potentially hazardous 
to EartrQ. One of these potentially hazardous asteroids, named 2012 DA14, was 
discovered in 2012. This asteroid will pass by Earth on February 15, 2013, coming 
closer to the Earth than satellites in geostationary orbij§. How close will 2012 DA14 



1 See http://neo.jpl.nasa.gov/stats/ 

2 See http://neo.jpl.nasa.gov/neo/groups.html 

3 See article about 2012 DA14 posted March 6, 2012 on MSNBC.com 
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pass by Earth? A mere 17000 miles (27000 kmfl In cosmic terms, this close shave 
of 2012 DA14 with Earth in 2013 is a near-collision motion. 



3 The A/-Body Problem 



To model collision and near-collision motions we make some simplifying assump- 
tions and use Newton's Inverse Square Law of Gravity. We assume that all the bodies 
are idealized as particles with zero volume (i.e., as points), that no particle is torn 
apart by tidal forces, that the mass of each particle never changes, and that besides 
Newton's Law of Gravity there are no other forces acting on the bodies. Under these 
assumptions we would think of Shoemaker-Levy 9 as not being torn apart by tidal 
forces, but as colliding with Jupiter as a whole. 



3.1 The Equations 

The particles modeling the bodies move in three-dimensional Euclidean space 
which we denote by R 3 . For a positive integer N > 2, suppose there are N parti- 
cles with positions G R 3 and masses mj > 0, j = l,...,N. The distance between 
two of the particles is denoted by 



which is the standard Euclidean distance between two points in R 3 . The Newtonian 
Af-Body Problem is the system of second-order nonlinear differential equations 



where ' = d/dt for a time variable t and G = 6.6732 x 10~ n m 2 /s 2 kg. By an appro- 
priate choice of units of the q ; , we will assume that G = 1 because we are investi- 
gating the qualitative or geometric, rather than the quantitative, behavior of collision 
and near-collision motions. 

By the standard existence, uniqueness, and extension theory in differential equa- 
tions (see ]9] for example), the initial value problem 



nk = |q/-q*l, j^k, 



Gmjm k (q k - qj) 



..,N, 



mjm k {q k -qj) 



<U(fo)=q?, q»=q?, 



(1) 



has a unique solution 



q(0 = (qi (0. — ,q*(0) 



4 See article about 2012 DA14 posted March 8, 2012 on Earthsky.org 
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defined on a maximal interval of definition (f~,f + ) as long as r,-j ^ for all j ^ k 
at t = tQ. Such a solution q(f ) describes a motion of the AT particles. 

Not every initial value problem (Q~|i will have a solution q(f ) with t~ = — °° and 
t + = oo. A solution with either r ~ > — °° or f + < °° experiences a singularity at the 
finite endpoint of its maximal interval of definition. The notion of a singularity is 
addressed in Section PTTl 



3.2 Integrals 

An integral of motion of the Newtonian A^-Body Problem is a differentiable function 
F of the position q and/or the velocity q' and/or the masses m = (mi,..., rn.fi) sucn 
that 

jF(q(t),q'(t),m) = 0, t€{t~,t + ). 

Along a solution q(f ) an integral F of motion satisfies 

F(q(t),q'(t),m)=F(q(t ),q'{t ),m), r <t <t + , 

i.e., the value of F is constant along the solution. The Newtonian A^-Body Problem 
has ten known integrals of motion. 

The translation invariance of the equations of the Newtonian A^-Body Problem 
give rise to 6 integrals of motion. With M = DyLi m j, three of these are given by the 
components of the center of mass vector 

^ m 

C = — V 111:0;, 

M p\ JHJ 

and three more are given by the components of the linear momentum vector 

1 N 
L = — V fw.-q';. 
M p x jHj 

Typically, both of these are set to so that the relative motion of the Af particles is 
emphasized. 

The rotational symmetry of the equations of the Newtonian A^-Body Problem 
give rise to 3 more integrals of motion. These integrals are given by the components 
of the angular momentum vector, 

N 

A = £m 7 q y xq}. 

7=1 

The angular momentum plays a key role in understanding collisions in the A^-Body 
Problem, as we will see later on. 
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There is one more integral of motion of the Newtonian A^-Body Problem. The 
self-potential (or negative of the potential energy) is 

„ mjm k 

j<k r ik 

The kinetic energy is 

K = \ t m J<-<- 

The total energy 

H=K-U 

is an integral of motion for the Newtonian A^-Body Problem. 

In the late 1800's, the mathematical strategy for "solving" the Newtonian N- 
Body Problem was to find enough "independent" integrals of motion |25l . This 
would implicitly give each solution as the curve of intersection of the hypersur- 
faces corresponding to the integrals of motion. Each solution q(f ) is a curve in R 6N . 
However, the intersection of the hypersurfaces of the ten integrals of motion gives a 
6A^ 10 > 1 dimension hypersurface in R 6 ^, which is not a curve! The ten known 
integrals of motion are independent of each other (one is not a function of the oth- 
ers), and are algebraic functions of positions, velocities, and masses. Are there any 
more algebraic integrals of motion? This was answered a long time ago in 1887- 
1888 by Bruns Q. 

Theorem 1. There are no algebraic integrals of motion independent of the ten 
known integrals of motion. 

Consequently, new integrals of motion, if any, can not be algebraic! In 1893, New- 
comb [20] lamented that no additional integrals had been found to enable the im- 
plicit solution of the 3-Body Problem. It is well-known that the Newtonian 2-Body 
Problem can be solved imi 3licitl>0, but all attempts to solve the A^-Body Problem 
with N>3 have been futilfl 

Typically then the solution q(f ) of the initial value problem (Q~|i is estimated nu- 
merically. From the constant total energy H along a solution q(f ), we observe that 
if any of the distances get close to 0, i.e., at least two of the particles are near 
collision, the self-potential becomes large and the kinetic energy becomes large too. 
The latter implies that the velocity of at least one of the particles becomes large, and 
the linear momentum L along q(f ) implies that the velocity of at least two particles 
becomes large. In particular, from the equations of the Newtonian Af-Body Problem, 
the particles that are near collision are the one with the large velocities. These large 
velocities presents problems for the numerical estimates of such a solution. 



5 See en.wikipedia.org/wiki/Gravitational_two-body_problem 

6 Karl Sundman did solve the 3-Body Problem when A ^ by convergent power series defined for 
all time, but the series converge too slowly to be of any theoretic or numerical use 1251 . 
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Rather than solving the A^-Body Problem for all of its solutions by finding enough 
independent integrals of motion, it is better to examine special solutions with par- 
ticular features. The simplest solutions to find are equilibrium solutions, where the 
position qj(t ) of each particle is constant for all time. But the Newtonian A^-Body 
Problem has none of these (see p. 29 in fl8l ). The next simplest solutions are peri- 
odic solutions, i.e., there exist T > such that q(f + T) = q(f ) for all t G R. These 
are part of the larger collection of solutions q(r ) with t ~ = — °° and t + = °° that are 
bounded. Such solutions must have a particular total energy (see p. 160 in |25l ). 

Theorem 2. If a solution q(f ) of the Newtonian N-Body Problem exists for all time 
and is bounded, then the total energy H < 0. 

Consequently, any periodic solution q(f ) of the Newtonian A^-Body Problem must 
has negative total energy. This is why in the search for periodic solutions the total 
energy is always assigned a negative value. 



3.4 Stability 

A periodic solution q(f ) of the Newtonian A^-Body Problem gives a predictable fu- 
ture: we know with certainty what the positions of the N particles will be at any 
time t > 0. But what if our measurements of the initial conditions q(0) and q'(0) 
are slightly off? A solution q(f ) with initial conditions near q(0) and q'(0) will stay 
close to q(f ) for a time, by a property of solutions of initial value problems called 
continuity of solutions with respect to initial conditions (see [9]). But if it stays close 
for all t > 0, we think of q(f ) as being stable. 

To quantify this notion of stability for a periodic solution, we use a Poincare 
section which is a hyperplane S containing the point (q(0),q'(0)) that is transverse 
to the curve (q(f ),q'(f))- If x = (q(0),q'(0)) is a point on S near the (q(0),q'(0)), 
then P(x) is next point where the the curve (q(f),q' (f)) intersects and P 2 (x) is 
the next point, and so on. The initial condition x° = (q(0),q'(0)) is a fixed point of 
this Poincare map P from S to S, i.e., P(x°) = x°. 

Definition 1. The periodic solution q(f ) is stable if for every real e > there exist a 
real 8 > such that \P k (x) -x°| < £ for all k = 1,2,3, ... , whenever |x - x°| < 8. 

When q(f ) is not stable, there are solutions which start nearby but eventually move 
away from q(f ), and we say that q(f ) is unstable. 

Showing directly that q(f ) is stable or unstable is very difficult. Instead, the re- 
lated concept of linearized stability is investigated, at least numerically. The deriva- 
tive of the Poincare map at the fixed point x° is a square matrix DP(x°). 



For an illustration of this see en.wikipedia.org/wiki/Poincare_map 
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Definition 2. A periodic solution q(f ) is 

1. spectrally stabl^if all the eigenvalues of DP(x°) have modulus one, and is 

2. linearly unstable if any eigenvalue of DP(x°) has modulus bigger than one. 

In 1907, Liapunov |fl5l established a connection between the stability of Definition 
Q]and the linearized stability of Definition [2] 

Theorem 3. If a periodic solution q(f) is stable then it is spectrally stable, and if 
q(t ) is linearly unstable, then it is unstable. 

If a periodic solution is shown numerically to be linearly unstable, then by Theorem 
3 the periodic solution is unstable. On the other hand, if a periodic solution is shown 
numerically to be spectrally stable, it may be stable or unstable. Examples exist with 
spectrally stable fixed points of maps like P that are unstable (see l28l ). 

The notion of stability for a non-periodic solution, such as the motion of the 
sun and planets in the solar system, is harder to grasp. Here is a sampling of the 
history and opinions on this stability problem. In 1891, Poincare commented that 
the stability of the solar system had at that time already preoccupied much time and 
attention of researchers (see p. 147 in iPTOl ). In 1971, Siegel and Moser lamented that 
a resolution of the stability problem for the TV-Body Problem would probably be in 
the distant future (see p. 219 in 11281 ). In 1978, Moser noted that the answer to the 
stability of the solar system was still not known (see p. 127 in ITOl ). In 2005, Saari 
stated that a still unresolved problem for the TV-Body Problem is that of stability 
(see p. 132 in fl25l ). Meyer, Hall, and Offin commented how little is known about 
the stability problem and how difficult it was to get (see p. 229 in |fT8l). 

In 1996, Diacu and Holmes suggested that the solar system should be consid- 
ered stable (in a weak sense) if no collisions occur among the sun and the planets, 
and no planet ever escape from the solar system (see. p. 129 in 1 10|). In this weak 
sense of stability, the solar system is stable for the next few billion years according 
to numerical work of Hayes fTTI in 2007. Much longer-term numerical studies of 
the solar system by Batygin and Laughlin [6] in 2008 using small changes in the 
initial conditions suggest that Mercury could fall into the sun in 1.261Gyo or that 
Mercury and Venus could collide in 862Myi0 and Mars could escape from the solar 
system in 822Myr. The Newtonian V-Body Problem thus suggests that in the near 
future, the Solar System should be free of collisions of planets and the Sun, with 
no planets escaping the solar system. But this still leaves open the possibility that 
smaller objects, such as asteroids and comets, could collide with any of the planets 
in the short and long term. Recall that there are nearly 9000 of those near-Earth as- 
teroids to consider, with 2012 DA14 making its near-collision approach with Earth 
on February 15 of 2013. 



There is a more restrictive notion of spectral stability known as linear stability that requires 
additional technical conditions on the square matrix DP(x°). 

9 Gyr means giga-year or 1 ,000,000,000 years 

10 Myr means mega-year or 1,000,000 years 
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4 Collisions 

Either in the short term of the long term, collisions put a wrench into the question 
of any notion of stability. Why should a solution or any nearby solution of the New- 
tonian A^-Body Problem be defined for all time? Remember that Shoemaker-Levy 9 

has t + < °°! 



4.1 Singularities 

Collisions are one of the two kinds of singularities in the Newtonian A^-Body Prob- 
lem. The solution q(f ) of initial value problem (Q]) is real analytic (i.e., a convergent 
power series) on an interval (?o — <5,fo + 8) for some 8 > 0, as long as ^ for 
all j ' ^ k at fo. By a process called analytic continuation (see for example |fT6j|), the 
interval (fo — 8 , fo + 8 ) can be extended to the maximal interval (f ~ , f + ) . 

Definition 3. A singularity of the Newtonian A^-Body Problem is a time f = f + or 
f ~ when f + < °° or f ~ > — °°. 

In 1897, Painleve Il22l characterized a singularity of the Newtonian A^-Body 
Problem, using the quantity 

rminO) =mirir,*( / ) 

#* 

determined by a solution q(f ). 

Theorem 4. A singularity for the Newtonian N-body Problem occurs at time t = t* 
if and only if V m j n (f) — > as t — > t *. 

An understanding what this means is obtained by considering the collision set 
4 = |J{q:q. = q fc }c(RT ) 

which is the set of points where two or more of the Af -particles occupy the same 
position. Painleve's characterization means that q(f) approaches the collison set, 
i.e., 

q(f) -> A as t ->t* 

when f* is a singularity of the Newtonian A^-Body Problem. Painleve's chararteriza- 
tion introduces two classes of singularities. 

Definition 4. A singularity f * of the Newtonian A^-Body Problem is a collision sin- 
gularity when q(t) approaches a specific point of A as f — > f *. Otherwise the singu- 
larity f * is a non-collision singularity. 
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Only collision singularities can occur in the Newtonian 2-Body Problem because 
it can be implicitly solved. In 1897, Painleve l22l showed that only one other New- 
tonian A^-Body problem has only collision singularities. 

Theorem 5. In the 3-Body Problem, all singularities are collision singularities. 

Unable to extend his result to more than 3 bodies, Painleve conjectured that there ex- 
ist non-collision singularities in the Newtonian 4 or larger Body Problem. In 1992, 
Xia IT341 mostly confirmed Painleve's conjecture, giving an example in the Newto- 
nian 5-Body Problem. 

Theorem 6. There exist non-collision singularities in the N-Body Problem for N > 
5. 

That leaves unresolved the question of the existence of non-collision singularities in 
the Newtonian 4-Body Problem. 

An understanding of what a non-collision singular looks like is obtained by con- 
sidering one-half of the polar moment of inertia of the Newtonian A^-Body Problem: 

1 N 

This scalar quantity measures the "diameter" of the N particles in the Newtonian 
A^-Body Problem. In 1908, von Zeipel 071 characterized a collision singularity in 
terms of the polar moment of inertia. 

Theorem 7. A singularity of the Newtonian N-Body Problem at t = t* is a collision 
if and only if I is bounded as t t*. 

This implies that for a non-collision singularity, at least one of the Af -particles has 
to achieve an infinite distance from the origin in just a finite time. This is a rather 
strange thing for Newton's Law of Gravity to predict. On the other hand, by The- 
orem [7] for a collision singularity, all of the positions of the N particles remain 
bounded at the moment of the singularity. 

A total collapse is an example of a collision singularity in the A^-Body Problem 
for which all N particles collide at the same point at the singularity t*. For a solution 
q(f ), the quantity 

»max = maxr / j t (?) 

characterizes a total collapse: a total collapse occurs at t * if and only if 

fmaxO) ->0 as f-M*. 

There is a relationship between total collapse and the angular momentum that was 
known by Weierstrass and established by Sundman (see l25lD . 

Theorem 8. If A ^ 0, then r max (f) is bounded away from zero. 

This does not preclude the collision of less than N particles when A ^ 0, as will be 
illustrated for certain Newtonian A^-Body Problems in Section[5] 
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4.2 Improbability 

Recall that there are 1306 potentially hazardous near-Earth asteroids. What are the 
chances that Earth will be hit by a near-Earth asteroid, or Jupiter will be hit by 
another comet? Well, it depends on the arrangement of the particles. 

Definition 5. A solution q(f ) is called collinear if the N particles always move on 
the same fixed line in R 3 . Otherwise it is called non-collinear. 

Every collinear solution has zero angular momentum because q,-(i ) is parallel with 
q'j(t) for all t G (f~,f + ). In 1971 and 1973, Saari l23ll24l established the probability 
of collisions. 

Theorem 9. The probability that a non-collinear solution q(f) will have a collision 
is zero. Every collinear solution q(f ) has a collision. 

With collision singularities being rare for a non-collinear A^-Body Problem, why 
bother to study them? Diacu and Holmes (see p. 84 and p. 103 in iflOl ) argue for 
the study of collision singularities because without such a study, a complete under- 
standing of the Newtonian A^-Body Problem could not be achieved. In particular, 
solutions near collision singularities could behave strangely, and the probability of 
a solution coming close to a collision singularity is positive and thus can not be 
ignored. Understanding then the collision singularities enables an understanding of 
the near-collision solutions. 



4.3 Regularization 

Regularization is one method by which we can get an understanding of a collision 
singularity. To regularize a collision means to extend the solution beyond the col- 
lision through an elastic bounce without loss or gain of total energy in such a way 
that all of the solutions nearby have continuity with respect to initial conditions, i.e., 
they look like the extended collision solution for a time (see p. 104 and p. 107 in 
fTOlD . Regularization is typically done by a Levi-Civita type change of the depen- 
dent variables, and a Sundman type change of the independent variable (see J8)), 
that together removes the collision singularity from the equations. We illustrate this 
regularization in the simplest of the A^-Body Problems. 

In the Collinear 2-Body Problem (or Col2BP for short), the positions of the two 
particles are the scalar quantities q\ and q%. If x — qi — q\ is the distance between 
the particle with mass m\ at q\ and the particle with mass m-i at qi > q\, then the 
Col2BP takes the form 

x" = 4-A *>0, (2) 

x 1 

and the total energy takes the form 
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OT1OT7 , mi mi 

H= - ( ' (x 1 ) 2 -^^. (3) 
2(mi +m2) x 

As x — >• the two particles approach collision, and the total energy implies that the 
two particles collide with an infinite velocity, 

(x'f -> 00. 

To regularize the binary collision (or total collapse) in this problem, define a new 
independent variable s and a new dependent variable w by 

ds 1 2 

— = -, vv =*, 

where the former is the Sundman type change of the independent variable, and the 
latter is the Levi-Civita type change of the dependent variable. If ' = d/ ds, the 
second-order equation (O becomes 

w 2 [2ww - 2w 2 + (mi + m 2 )] = 0, (4) 

and the total energy (01 becomes 

2 2mim 2 . 2 

Hw = w —m\ni2- (5) 

ni\ +ni2 

As w — > 0, the second-order equation makes sense (no dividing by zero), and the 
total energy (0 implies that 

2 m l +m 2 

(w) 

which is a finite nonzero velocity! The collision singularity has been regularized. 

The regularized nonlinear second-order equation can actually be solved! 
Solving the total energy <(5j for 2(vv) 2 and substituting this into the second-order 
equation © gives 

(mi +nn)H 

w w 

2mini2 



2v? 



0. (6) 



This makes sense when w = 0, i.e., the moment of collision! For negative H, the 
linear second-order equatior{3 inside the square brackets in (|6| solves to give a 
real analytic stable periodic solution w(s) which experiences a collision every half 
period in terms of the regularized time variable s. The corresponding solution x(t) 
is periodic and experiences a collision once a period in terms of the original time 
variable t . This doubling of the number of collisions per period is because the change 
of dependent variable w 2 = x has w(s) "doubling" x(t) in that w(s) passes through 
twice a period, going from positive to negative and then negative to positive, while 
x(t) is positive except at collision where it is zero. 



11 This is a simple harmonic oscillator for H < whose solutions are in terms of cosine and sine. 
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The binary collision singularity in the Newtonian 2-Body Problem can be reg- 
ularized in a similar but more complicated way than what was done above for the 
Col2BP (see fl25l ). By Theorem[8] a solution of the 2-Body Problem with nonzero 
angular momentum does not experience a collision or total collapse. A nonzero 
angular momentum near-collision solution looks like the zero angular momentum 
collision solution^! The regularized 2-Body Problem provides good numerical es- 
timates of the motion because there are no infinite velocities! 



4.4 McGehee 

What about regularization of a triple collision, when three of the particles meet? In 
1974, McGehee ifTTl showed that regularization of a triple collision is in general 
not possibl^H- Starting close together, two solutions that approach a near triple col- 
lision can describe radically different motions after the near triple collision. This 
kind of behavior is known as "sensitive dependence on initial conditions," and is 
an antithesis of stability. Triple collisions present a numerical nightmare! By exten- 
sion, collisions with four or more particles present the same nightmare! So the only 
regularizable collisions are those that are essentially a binary collision. 



5 Results 

Spectrally stable periodic solutions have been found in Newtonian A^-Body Prob- 
lems with regularizable collisions for N > 3. Three of these situations discussed 
here are the Collinear 3-Body Problem (or Col3BP), the Collinear Symmetric 4- 
Body Problem (or ColS4BP), and the Planar Pairwise Symmetric 4-Body Prob- 
lem (or PPS4BP). There are other Newtonian A^-Body Problems where periodic 
solutions with regularizable collisions whose existence has been given analytically 
Il35l[36ll27l . some of whose stability (in the sense of Definition[T]i and linear stability 
(as defined in Definition [2]i has been numerically determined 151 |33] l35l [36l . 



5.1 Col3BP 

As a subproblem of the Newtonian 3-Body Problem, the Col3BP requires that the 
three particles always lie on the same line through the origin. The positions of the 



Binary Star Systems are known to exist in the Universe. The Newtonian 2-Body Problem pre- 
dicts stability for a Binary Star System, a collision-free solution that is bounded for all time. 
13 This is achieved by "blowing-up" the triple collision singularity and slowing down the motion 
as the particles approach a triple collision. This setting does allow for good numerical estimates of 
near-triple collisions. 
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three particles in the Col3BP are the scalars q\, qj_, and qj, which can be assumed to 
satisfy 

q\<qi< q-i- 

By Theorem |9j collisions always occur in the Col3BP. Because the three particles 
are collinear for all time, their angular is zero, and by Theorem [8] a total collapse 
is possibkQ in the CoOBP. In 1974, S.J. Aareth and Zare fT] showed that any two 
of the three possible binary collisions in the 3-Body Problem are regularizablj^l. 
In 1993, Hietarinta and Mikkola IT3l used Aarseth and Zare's regularization |T) to 
regularize the binary collisions q\ = qi and qi = q?, in the Col3BP. 

In 1956, Schubart fl26l numerically found a periodic orbit in the equal mass 
CoOBP of negative total energy in which the inner particle oscillates between bi- 
nary collisions with the outer particles. In 1977, Henon lfl2l numerically extended 
Schubart's periodic solution to arbitrary masses and investigated their linear stabil- 
ity. In 1993, Hietarinta and Mikkola (131 also numerically investigated the linear sta- 
bility of Schubart's periodic solution for arbitrary masses. Together they showed that 
Schubart's periodic solution is spectrally stable for certain masses, and linearly un- 
stable for the remaining masses. Hietarinta and Mikkola fl3\ further numerically in- 
vestigated the Poincare section for Schubart's periodic solution for arbitrary masses, 
showing when there is stability as described in Definition!]] In 2008, Moeckel |fl9l 
and Venturelli |32l separately proved the analytic existence of Schubart's solution 
when m\ = to 3 and TO2 is arbitrary. Only recently, in 201 1, did Shibayama l27l ana- 
lytically prove the existence of Schubart's periodic solution for arbitrary masses in 
the Col3BP. 

Schubart's period solution for the CoOBP is also a periodic solution of the 3- 
Body Problem, where in the latter the continuity with respect to initial conditions 
can be see for near-collision solutions. For example, Schubart's periodic solution 
for the nearly equal masses 

mi = 0.333333, m 2 = 0.333334, m 3 = 0.333333 

is spectrally stable. Considered in 3-Body Problem, Schubart's periodic solution for 
these mass values remains spectrally stable fT2l , and numerically the near-collision 
solutions in the Newtonian 3-Body Problem behave like Schubart's periodic solu- 
tion. It is therefore possible that in the 3-Body Problem, there are solutions near 
Schubart's periodic solution that are free of collisions and bounded for all time. 
Imagine, as did Henon 1121 . of Newton's Law of Gravity predicting a triple star 
system that is free of collisions and bounded for all time! 



Initial conditions leading to total collapse in the equal mass Col3BP are easy to realize: set 
q\ = — 1, Cj2 = 0, and qj, = 1 with the initial velocity of each particle set to 0. 
15 A good numerical model for the Sun-Jupiter-Shoemaker-Levy 9 or Earth-Moon-2012DA14 sit- 
uation is regularized 3-Body Problem of Aarseth and Zare. 
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5.2 ColS4BP 

As a subproblem of the Newtonian 4-Body Problem, the ColS4BP requires that the 
four particles always lie on the same line through the origin. The positions of the 
four particles are the scalars q\, q 2 , q^, and q\ that satisfy 

<?4 = -qi, q3 = -q%, q\ >o, q 2 > o, 

and 

-qi < -qi <0<q 2 <qi 

with masses 

m l = 1 , 7772 = m > 0, '«3 = m, OT4 = 1 . 

The angular momentum for all solutions of the ColS4BP is zero because of the 
collinearity, and so by Theorem[8]a total collapse is possible. There are two kinds of 
non-total collapse collisions in the ColS4BP: the binary collision of the inner pair 
of particles of mass m each, i..e, q 2 = 0, and the simultaneous binary collision of the 
two outer pairs of particles, i.e., q\ — q 2 > 0. In 2002 and 2006, Sweatman ||30ll3Tl 
showed, by adapting the regularization of Aarseth and Zare JT|, that these non-total 
collapse collisions in the ColS4BP are regularizable. 

Sweatman [30, 31] numerically found a Schubart-like periodic solution in the 
ColS4BP with negative total energy for arbitrary m where the outer pairs collide in 
a simultaneous binary collision at one moment and then the inner pair collides at an- 
other moment. He determined numerically that this Schubart-like periodic solution 
is spectrally stable when 

< m < 2.83 and m > 35.4, 

and is otherwise linearly unstable. In 2010, Bakker et al [{2) verified Sweatman's 
linear stability for the Schubart-like periodic solution using a different technique. In 
2011-2012, Ouyang and Yan |21|, Shibayama [[27] , and Huang ITT41 proved sepa- 
rately the analytic existence of the Schubart-like periodic solution in the ColS4BP. 



5.3 PPS4BP 

The PPS4BP has two particles of mass 1 located at the planar locations 

qi andq 3 = -q l5 
and two particles of mass < m < 1 located at the planar locations 

q 2 and q 4 = -q 2 . 
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The four particles in the PPS4BP need not be collinear, so that the angular mo- 
mentum need not be 0. Unlike the ColS4BP, total collapse can be avoided in the 
PPS4BP by Theorem[8]when the angular momentum is not zero. Like the ColS4BP, 
there are two kinds of non-total collapse collisions in the PPS4BP: simultaneous 
binary collisions when qi = q2 and q3 = q4, or when qi = q4 and q2 = q3; and 
binary collisions when qi = or when q2 = 0. In 2010, Sivasankaran, Steves, and 
Sweatman [29] showed that these non-total collapse collisions in the PPS4BP are 
regularizable. 

The Schubart-like periodic solution in the ColS4BP is also a periodic solution of 
the PPS4BP, where in the latter the continuity with respect to initial conditions can 
be observed for near-collision solutions. However, as shown by Sweatman PP . in 
the PPS4BP the Schubart-like periodic solution of the ColS4BP becomes linearly 
unstable for 

< m < 0.406, and 0.569 < m < 1.02 
as well as 2.83 < m < 35.4, while it remains spectrally stable for 

0.407 < m < 0.567 and m > 35.4. 

By long-term numerical integrations for the Schubart-like periodic solution as a so- 
lution of the PPS4BP, Sweatman 11311 showed that stability in the sense of Definition 
Q]is possible when 0.407 < m < 0.567 and when m > 35.4. It is therefore possible 
for these values of m that near Schubart's periodic solution there are collision-free 
solutions of the PPS4BP that are bounded for all time. 

In 2011, adapting the regularization of Aarseth and Zare HI to simultaneous 
binary collisions, Bakker, Ouyang, Yan, and Simmons y] proved the analytic exis- 
tence of a non-collinear periodic solution in the equal mass PPS4BP. This periodic 
solution has zero angular momentum, negative total energy, and alternates between 
a simultaneous binary collision of the symmetric pairs in the first and third quadrant 
where qi = q2 and q3 = q4, and the simultaneous binary collision of the symmet- 
ric pairs in the second and fourth quadrants where qi = q4 and q2 = q3- Bakker, 
Ouyang, Yan and Simmons O then numerically extended this non-collinear peri- 
odic simultaneous binary collision solution to unequal masses < m < 1. In 2012, 
Bakker, Mancuso, and Simmons |4} have numerically determined that the non- 
collinear periodic simultaneous binary collision solution is spectrally stable when 

0.199 < m < 0.264 and 0.538 < m < 1 

and is linearly unstable for the remaining values of m. Long-term numerical inte- 
grations of the regularized equations done by Bakker, Ouyang, Yan, and Simmons 
(3) suggest instability when 0.199 < m < 0.264 and stability when 0.538 < m < 1 
in the sense of DefinitionQ] For these latter values of m could the near-collision so- 
lutions in the PPS4BP that look like the non-collinear periodic simultaneous binary 
collision solution be collision-free and bounded for all time? 
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Both the ColS4BP and the PPS4BP are subproblems of the Newtonian 4-Body 
Problem, where the non-total collapse collisions in the former two problems are 
regularizable. What is not known is how to, if possible, regularize binary collisions 
and simultaneous binary collisions in the Newtonian 4-Body Problem within one 
coordinate system If such a regularization is possible, then all of the periodic 
solutions thus known in the ColS4BP and PPS4BP would also be periodic solutions 
of the Newtonian 4-Body Problem and the investigation of their stability and linear 
stability in the Newtonian 4-Body Problem could begin. With more possible per- 
turbations of initial conditions in the Newtonian 4-Body Problem as compared with 
the PPS4BP, a loss of spectral stability could indeed happen as it did with going 
from the ColS4BP to the PPS4BP. But some of the spectral stability might survive 
passage from the PPS4BP to the Newtonian 4-Body Problem, giving the possibility 
of near-collision solutions that are collision-free and bounded for all time. 
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